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Abstract. Supermassive black hole binaries (SMBHBs) with masses in the mass 
range ~ (10^-10^) Mq/(1 + z), produced in galaxy mergers, are thought to 
complete their coalescence due to the emission of gravitational waves (GWs). 
The anticipated detection of the GWs by the future Laser Interferometric 
Space Antenna (LISA) will constitute a milestone for fundamental physics and 
astrophysics. While the GW signatures themselves will provide a treasure 
trove of information, if the source can be securely identified in electromagnetic 
(EM) bands, this would open up entirely new scientific opportunities, to probe 
fundamental physics, astrophysics, and cosmology. We discuss several ideas, 
involving wide-field telescopes, that may be useful in locating electromagnetic 
counterparts to SMBHBs detected by LISA. In particular, the binary may produce 
a variable electromagnetic flux, such as a roughly periodic signal due to the orbital 
motion prior to coalescence, or a prompt transient signal caused by shocks in the 
circumbinary disk when the SMBHB recoils and "shakes" the disk. We discuss 
whether these time- variable EM signatures may be detectable, and how they can 
help in identifying a unique counterpart within the localization errors provided 
by LISA. We also discuss a possibility of identifying a population of coalescing 
SMBHBs statistically, in a deep optical survey for periodically variable sources, 
before LISA detects the GWs directly. The discovery of such sources would 
confirm that gas is present in the vicinity and is being perturbed by the SMBHB 
— serving as a proof of concept for eventually finding actual LISA counterparts. 



1. Introduction 

The anticipated detection by LISA of gravitational waves emitted during the 
coalescence of supermassive black holes (SMBHs) in the mass range ^ (10'^- 
lO'') M0/(1 -I- z) will constitute a milestone for fundamental physics and astrophysics. 
While the GW signatures themselves are a rich source of information, if the GW source 
produces electromagnetic (EM) radiation, and if the object can be securely identified in 
EM bands, this would open up entirely new scientific opportunities. The simultaneous 
study of photons and gravitons from a single source could probe fundamental aspects 
of gravitational physics [8, 15]. The GW sources can also be used as self-calibrated 
standard sirens [26], and cosmological parameters can be determined if the source 
redshift is identified [11, 13]. Finally, for many events in the above mass and redshift 
range, LISA will be able to measure the masses and spin vectors of the SMBHs, their 
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orbital parameters, and their luminosity distance, to a precision unprecedented in any 
other type of astronomical observation [17, 28]. If a counterpart is known, then the 
Eddington ratios and other attributes of black hole accretion physics can be studied 
in exquisite detail [13, 15]. 

Motivated by the above possibilities, several recent studies have addressed the 
question of whether finding a counterpart will be feasible, given LISA^s localization 
errors [11, 13, 14, 18, 15]. The crucial uncertainty, of course, is the nature (luminosity, 
spectrum, and time-evolution) of any EM emission produced by coalescing SMBHBs 
during the GW-inspiral stage. Such emission would have to be related to gas in the 
vicinity, and possibly accreting onto the coalescing SMBHs. The gas around the BHs 
would likely settle into a rotationally supported, circumbinary disk. In a geometrically 
thin disk, the torques from the binary create a central cavity, nearly devoid of gas, 
within a region about twice the orbital separation [2] (for a nearly equal-mass binary), 
or a narrower gap around the orbit of the lower-mass BH in the case of unequal masses 
q = M1/M2 <C 1 and larger orbital separations [1]. In the latter case, the lower-mass 
hole could "usher" the gas inward as its orbit decays, producing a prompt and luminous 
signal during coalescence. In the former case, residual gas flow into the cavity and 
onto the BHs, such as suggested in numerical simulations [3, 20, 7], may still produce 
non-negligible EM emission. Around the time of coalescence, the gravitational waves 
shear the circumbinary gas and could brighten its emission dctectably [16]. Finally, 
SMBHBs recoil at the time of their coalescence due to the emission of GWs, at speeds 
up to 4,000 km s~^ (e.g. ref. [5] and references therein). The gas disk will respond 
promptly (on the local orbital timescale) to such a kick, which may produce prompt 
shocks, and transient EM emission after coalescence [19]. J 

The complex processes involved in ultimately producing any EM emission remain 
poorly understood, and the level of the luminosity produced, as well as its spectrum 
and time-evolution, are essentially unknown. However, any emission during and 
promptly following the inspiral stage is likely to be variable. In this contribution, we 
discuss three issues related to using variability to identify EM counterparts. Can 
LISA events be localized to within the field of view of astronomical instruments 
(several deg^), hours to weeks prior to coalescence (§ 2)? What is the response of 
the circumbinary gas to the gravitational recoil ( "kick" ) of the SMBHB at coalescence 
(§ 3)? Can we identify coalescing SMBHBs before the launch of LISA, as variable 
sources, due to periodic perturbations in the circumbinary gas (§ 4)? 

2. Monitoring the 3D LISA Error Box Prior To Coalescence 

The first and most fundamental question in searching for any EM counterpart by 
looking for variable emission, during the final stages of coalescence, is the accuracy 
to which the LISA source can be localized at various look-back times prior to the 
coalescence. (Coalescence is taken to occur when the binary separation reaches the 
innermost stable circular orbit; ISCO). Here we present time-dependent localization 
errors, obtained by the Harmonic Mode Decomposition method [14]. This technique 
uses the restricted, post Newtonian approximation for the GW waveform, and applies 
the Fisher matrix technique to the Fourier transform of the waveform, to forecast 
parameter uncertainties. Orbits are assumed to be circular, and spins are neglected 

X Such brightening due to kick-induced disk heating can persist on much longer timescales, and could 
produce detectable emission - leading to recent proposals that this can help identify a population of 
such sources before EISA's launch. This possibility will be briefly discussed further in § 4 below. 
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Figure 1. Evolution witli pre-ISCO look— back time, t{, of the expected LISA 
source localization uncertainties. Left Panel: Sky position errors (major axis 2ajv, 
minor axis 26 at, and equivalent diameter, 2rfj = 2\/ajv5jVi of the error ellipsoid). 
An equal-mass binary is shown as an example, with Mi = M2 = IO^Mq at 2 = 1. 
Best, typical, and worst cases (among sources with random orientation) represent 
the 10%, 50%, and 90% levels of cumulative error distributions. The absolute 
errors are typically small enough to fit in the FOV of a wide— field instrument in 
the last few weeks, allowing a world-wide monitoring campaign. The horizontal 
line shows a diameter of 3.57°, which corresponds to localizing the source to within 
10 deg^. Right Panel: Luminosity distance di^{z) errors for the same binary (10%, 
50%, and 90% levels of cumulative error distributions as in the left panel). The 
horizontal line delineates the level of weak lensing uncertainties (adopted from 
ref. [13]). Knowing the luminosity distance to the accuracy of a few percent will 
allow the counterpart search to focus on candidates in a narrow redshift slice. 

(this is justified until the last day or so of the merger; see [18]). The 17-dimensional 
parameter space describing the general binary inspiral is split into "slow" and "fast" 
parameters, based on the time-scales on which they modulate the waveform, and the 
two sets are assumed to be decoupled. The angular coordinates of the source, and its 
luminosity distance - representing the 3D localization of the source - are "slow" Fisher 
parameters, and information on these parameters are derived from the annual motion 
of LISA around the Sun. The reader is referred to [14] for a full list of assumptions 
and details about the method. 

The time-evolving Icr errors on the two-dimensional sky position are shown for an 
equal-mass binary, Mi = M2 — 10^ M©, at redshift z = 1 in Figure 1 (corresponding 
roughly to the optimal choice of mass/redshift combination; other masses and redshifts 
yield poorer localization) . The HMD method, by construction, approximates the sky 
position errors by ellipses. Figure 1 shows the gradual improvement in localization, 
in the form of the major axis (20^?), minor axis (2b]\;) and equivalent diameter 
{2rN = -\/4a7v&Ar). Figure 1 displays results for three separate cumulative probability 
distribution levels, 90%, 50%, 10%, so that 10% refers to the best 10% of all events, 
as sampled by the random distribution of five angular parameters. The evolution 
of errors scales steeply with look-back time for t{ > 40 days. For smaller look-back 
times, errors essentially stop improving in the "worst" (90% level) case, improve with 
a relatively shallow slope for the "typical" (50% level) case, and improve more steeply 
in the "best" case (10% level among the realizations of fiducial angular parameters). 
Similar evolutionary trends are seen for the luminosity distance ^^(z) errors (shown 
in the right panel in Fig. 1). 

Figure 2 displays contours of fixed "advance warning time" for typical (50%) 
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Figure 2. Contours of advance warning times in the total SMBHB mass (M) vs. 
redshift (z) plane, for typical events with SMBH mass ratio q = M1/M2 = 1. The 
contours trace the look-back times at which the major axis (2a) of the localization 
error ellipsoid first falls below 3.57°. The contours are logarithmically spaced in 
days, and 10 days is highlighted with a thick curve. 

events, adopting a 10 deg^ FOV as a reference. The contours are logarithmically 
spaced, with solid contours every decade and the shaded region highlights the (Af, z) 
region where at least 10~day advance notice will be available. This figure shows that 
10 day advance warning to cover the full error ellipsoid with a single LSST pointing is 
possible for a range of masses and source redshifts, up to Af ^ 3 x IO^Mq and z ~ 1.7. 

The results shown in Figures 1 and 2 suggest that it will be possible to identify, 
prior to merger, a small enough region in the sky where any prompt electromagnetic 
(EM) counterpart to a LISA inspiral event will be located. Given sufficient advance 
notice, it will then be possible to trigger a world-wide monitoring campaign, to search 
for EM counterparts as the merger proceeds (and also during the most energetic 
coalescence phase). 

The several square-degree field will, of course, contain a very large number of 
sources (a few x 10^ galaxies in total, at the limiting optical magnitude of ~ 27 mag 
that may be relevant; e.g. ref. [21] and discussion below). Having predictions for the 
spectrum and time-dependence of a coalescing SMBH binary, and therefore knowing 
what to look for, would obviously greatly help in identifying the counterpart, possibly 
allowing an identification even post-merger. Such a prediction, however, would require 
understanding the complex hydrodynamics and radiative properties of circumbinary 
gas at the relevant small separations of a few to a few thousand Schwarzschild radii. 
This is a notoriously difficult problem even in the much simpler case of steady accretion 
onto a single SMBH [22]. This suggests that the best strategy may be an "open- 
minded" search for any variable signatures prior and during coalescence. 

There will be several ways, however, to cut down on the list of possible 
counterparts, using the photometric redshifts of the candidates (restricting the search 
to within a narrow, Sz <J few percent, redshift slice; see the right panel in Fig. 1), the 
expected luminosity of the source, and the other parameters of the binary provided by 
LISA (for example, a variable EM signal may be much more likely if the spin and the 
orbital angular momenta are known to be aligned, as this may indicate the presence 
of circumbinary gas). These, and several other possible cuts are discussed further in 
ref. [15]. 
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3. Prompt Shocks in the Gas Disk Around a Recoiling Binary 

The recent break-through in numerical relativity has allowed a direct computation of 
the linear momentum flux ( "kick" ) produced during the coalescence of a SMBH binary. 
Such kicks may help produce prompt EM counterparts to GW sources detected by 
LISA. If the SMBHB is surrounded by a circumbinary gas disk, the disk will indeed 
respond promptly (on the local orbital timescale) to such a kick. If this results in 
warps or shocks, the disturbed disk could produce a transient EM signature [23]. 
As discussed above, the final sky localization uncertainty from LISA is typically a 
few tenths of a square degree, containing a large number of sources; monitoring this 
area for transient events after the merger may be another method to securely identify 
counterparts. 

We investigated the response of a circumbinary disk to the kick [19]. We 

adopted the following simplified picture for the disk around a fiducial, equal-mass. 
M1 + M2 = 10^ Mq binary. The disk has an inner edge at lOOrs (Schwarzschild radii), 
inside which it is empty (with the gas evacuated due to torques from the binary; e.g. 
ref. [23]) and an outer edge at 10,000rs. Outside of this radius, there may still be 
gas, but it may be unstable to fragmentation, and it will, in any case, evolve slowly 
(the behavior of this gas will then not be relevant for a LISA counterpart search, but 
if emission is produced in this gas, it could help identify SMBHBs independently; see 
discussion in the next section). The scale-height and temperature at the inner edge is 
h/r = 0.46 and T = 1.7 x lO^K, respectively. The scale-height remains constant with 
radius out to 2,000rs, beyond which it increases nearly linearly {h cx 7-21/20^^ -jj^g 
temperature varies with radius as T cx r"^/^" (see ref. [19] for further details). 

The important features of such a disk (as well as other proposed variants of thin 
a-disks) are the following: (i) orbital motions in thin disks are supersonic, so that the 
gas is susceptible to shocks if disturbed; (ii) at the relevant radii outside lOOrg the 
viscous time-scale is long, and the orbits are near Keplerian; (iii) gas near the inner 
edge of disk is tightly bound to the kicked BHB (worbit ~ 2.1 x 10"* km s~^), but the 
outer edge (worbit ^ 2.1 x 10'^ km s~^) can be marginally bound, or even unbound, 
for large kicks (a rough condition for being bound is Worbit ^ Wkick), and (iv) the total 
disk mass within 10, OOOrg is much less than the BHB mass, which justifies ignoring 
the inertia of the gas bound to the BHB. 

For a quantitative assessment of the disk's response to the kick, we employed the 
following approximation: the disk particles are assumed to be massless, coUisionless, 
and initially on co-planar, circular orbits. The kick simply adds the velocity Wkick 
to the instantaneous orbital velocity of each particle (in the inertial frame centered 
on the BHB). We used N — 10^ particles, distributed randomly and uniformly 
along the two-dimensional surface of the disk. The kick velocity was varied between 
500 km s^^ < i^kick < 4,000 km s~^, and directed either perpendicular or parallel to 
the initial disk plane. Note that in both the perpendicular and parallel case, we start 
with a two dimensional particle distribution (i.e. an infinitely thin disk), but in the 
perpendicular case, we then follow the orbits in 3D. 

Figure 3 shows, as an example, a face-on view of the surface density of the disk 
90 days after a kick with iJ^ick = 500 km s~^ in the plane of the disk (left panel). 
The sharp, tightly wound spiral features clearly seen in the figure trace the locus of 
points where particles cross each other, corresponding formally to a density caustic. 
The spiral caustic first forms at ^ 30 days, and then propagates outward at a speed 
of « 500 km s~^, so that the outermost caustic at time tc is located at Tc ^ tci'kick 




Figure 3. Top Panel: The top view of tlie surface density of a disk around 

an Ml + M2 = 10^ Mq BH binary, which recoiled within the disk plane at the 
velocity of tJ^ick = 500 km s^^, oriented vertically upward in the diagram. The 
disk is initially assumed to have an inner edge at r^^ = lOOrs (Schwarzschild 
radii), and is shown here out to a radius of Hn = 5,000rs, at a time t = 90 
days after the kick. A tightly wound, outward-moving spiral caustic develops 
in about a month. The dark/light shades correspond to regions of high/low 
density, the low-density regions being similar to the initial surface density, and 
the high-density regions about 10 times overdense. Bottom Panel: Aerial view 
of a disk with a kicked SMBHB, as in the left panel, except here the kick is 
oriented perpendicular to the disk, and the snapshot is taken at i = 1 week. 
For visual clarity, the graphic contrast was increased relative to the left panel, 
and the diagram stretched by a factor of 10 along the axis perpendicular to the 
disk plane. The density distribution is azimuthally symmetric, and while there 
are mild concentric density fluctuations, strong enhancements (i.e. caustics) were 
found to develop only after a delay of Ri one year. 

(this behavior can be roughly understood using epicycles; refs. [19, 25]). The right 
panel in Figure 3 shows, for comparison, the aerial view of the 3D particle density one 
week after a kick with the same velocity, but perpendicular to the disk. The density 
profile in this case remains azimuthally symmetric, but still develops concentric rings 
of density fluctuations (although we find the density enhancements are much weaker, 
at the ten percent level). 

These results suggest that strong density enhancements can form promptly after 
a supersonic kick in the plane of the circumbinary disk, within a few weeks of the 
coalescence of a ~ 10^ M© BHB. Because the disk is cold, and caustics are formed 
when particles first cross each other along their orbits, this implies that corresponding 
shocks could occur in a gas disk. For hydro dynamical shocks to occur within a 
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finite-pressure gas, the relative motions Vc between the neighboring particles that 
produce the caustic must exceed the sound speed. At the outermost radius where 
the disk is marginally bound to the BHB, one expects Vc ~ t^kick ^ I'orbHi relative 
motions will be slower further inside. The relative speed should roughly correspond 
to covering the cpicyclic amplitude ~ (I'kick/'f orbit )''c in the caustic-formation time 
tc ~ rc/vkick, yielding Vc ~ i^kick/^orbit- For i^kick = 500 km s~\ this predicts 
Vc ~ 25 km s~^(r/1000rs)^/^; we have verified in our simulations that particles cross 
the caustics at speeds within ~ 30% of this predicted value. Compared with the sound 
speed Cg « 25 km s~^(r/1000rs)~^/^°, this suggests that the density waves produced 
by the kick in the gas beyond 700rs will indeed steepen into shocks. We also found 
that the inclination of the kick may be important in determining the strength and 
timing of such shocks. 

The nature of the emission resulting from the shocks or density enhancements will 
have to be addressed in future work, by incorporating the effects of hydrodynamics, 
computing the heating rate at the location of the spiral shocks, and modeling the 
overall disk structure and vertical radiation transport. § However, in the fiducial 
case discussed above, if wc assume that the shocked gas is heated to temperatures 
corresponding to Vc, and tc is the time-scale on which the corresponding thermal 
energy is converted to photons, then we find that the luminosity may be a small but 
non-negligible fraction, 0.2 percent to 5 percent of the Eddington luminosity of the 
central BH, and would increase with time roughly as Lkick cx: t'^. This suggests that the 
afterglows may be detectable, at least for nearby BHs and/or for the most massive BHs 
in LISA'S range (which extends up to Mbh ~ lO^Mo). Wc may also speculate on the 
spectral evolution of the "kick after-glow" , based on the characteristic shock velocity 
at each radius. We find that the luminosity is dominated by the outermost shocked 
shells, with the spectrum peaking at the characteristic photon energy corresponding to 

^-^shock 

oc oc v^^l^^ DC r. The shocks could therefore result in an afterglow, starting 
from 700rs/?;kick ~ 50 days, first peaking in the UV band (~ 3eV or ~ O.S/um, 
corresponding to ^ 25 kms^-'^), and then hardening to the EUV/soft X ray 50cV) 
range after ~two years. The detection of such an afterglow would help identify EM 
counterparts to LISA events. 

4. Secirching for Gravitational Binary Inspirals Before LISA 

Numerical simulations suggest that the central cavity of the circumbinary disk is not 
completely empty, and that there can be non-negligible gas inflow into this cavity from 
the outside disk [3, 20, 7]. The perturbations of the circumbinary gas by the rotating 
quadrupole potential of the binary, in fact, appears to impose large fluctuations on 
this inflow rate, tracking the orbital period [20, 7]. We must emphasize that the level 
and nature of EM emission, associated with this inflow (or with other effects from the 
gas in the vicinity of the binary during the late stages of binary evolution), remains 
essentially unknown. Furthermore, even in the most optimistic scenario, in which a 
clearly periodic emission is indeed produced, at a significant and detectable flux level, 
there remains the (potentially severe) observational challenge of distinguishing such 
variability from other possible sources of periodic variations. 

Despite these uncertainties and caveats, however, it is interesting to ask the 

§ Our preliminary work, based on hydrodynamical simulations, does indicate that gas disks with 
realistic pressure profiles still develop strong shocks [6] . 
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following question: if, indeed, non-negligible emission is produced during the late 

stages of binary evolution, and the luminosity varies periodically on the orbital time- 
scale, could such periodic sources possibly be identified in EM surveys, even before 
LISA becomes operational? 

As the orbit of each individual binary shrinks, its orbital period, and therefore, by 
assumption, its variability timescale decreases. The evolution of a binary embedded 
in a thin circumbinary disk generically proceeds through three distinct stages (see 
ref. [10]). (i) First, the binary is strongly coupled to the circumbinary disk and is 
driven by viscosity (analogous to "disk-dominated" planetary migration), and the 
radius of the gap follows the binary. As the binary separation shrinks below ~ lO^rs 
(lO^rs) for mass ratios q = M1/M2 ^ 1 {q ^ 0.01), the local disk density starts 
to dominate over the binary mass, and the binary evolves more slowly, according 
to so-called "secondary-dominated" Type-II migration. During this stage, the GW 
emission is negligible, (ii) Later, within the radius ^ 500rs, the binary starts to be 
driven primarily by GWs but the radius of the gap can still follow the binary. (Hi) 
Finally, within ~ lOOrs the binary is entirely driven by GWs and the binary falls 
in much more quickly than the outer edge of the gap is able to move inward. The 
ordering of these events is valid for a very broad range of binary and disk parameters. 
Note that the outcome of the gas inside the binary's orbit is left unspecified above 
(see rcf. [1] for a possible outcome for q <C 1). 

In the last, GW-driven regime, the binary spends a characteristic time Tqw at 
each orbital separation rorb < ''GW that scales with the corresponding orbital time 
as T oc ^Q^b- T^^^ incidence rate of sources that have similar inferred BH masses, 
and show near-periodic variability on the time-scale tvar ~ iorb, would then follow 
N oc tvLv- On the other hand, at larger separations, where the evolution is due to 

viscous processes, the incidence rate has a much flatter dependence. For near-equal 

25/51 7/12 

mass binaries, the scaling is between N oc tvar and cx tvar , depending on whether 
the opacity is dominated free-free absorption or electron scattering, respectively. For 
g <i; 1, the scaling is slightly steeper, between N oc t^l^ and N oc t^'^l^^ for free-free 
absorption or electron scattering opacity, respectively (see the Appendix in ref. [10]). 
In the left panel of Figure 4, we show that the time spent at each orbital period can 
be interestingly long, i.e. a non-negligible fraction of the expected lifetime of quasar 
activity of a ^ few x 10^ years. 

Can we detect the flux variations from these sources? Luminosity variations 
corresponding to a fraction /Edd ^ 0.01 of the Eddington luminosity (or, e.g., a 
periodic component with ampHtude i w 26 + 2.5 log[(/Edd/0.01)(AfBH/3 x lO'^M©)"^] 
magnitudes in the optical for BHBs at z = 2) would have not been found in existing 
variability surveys. However, as shown in Figure 4, a dedicated, deep survey of a 
^deg^ area, possibly with existing instruments, could detect a population of ^several 
to ^several thousand such "transition sources" with a range of periods between 20 
weeks ^ t^ar ^ 60 weeks. In the right panel of Figure 4, we use the luminosity function 
of optical quasars, and assume that each optical quasar results from a SMBH SMBH 
merger, and, as the binary orbit decays, it goes through the evolutionary stages 
discussed above (see ref. [10] for more details). The figure then shows the required 
depth and area coverage for a multi year survey to detect such periodic sources, 
covering a factor of three range in period. The dependence in the abundance of these 
sources on their period - in particular, a switch from a flat, viscosity-driven power-law 
between N oc tiir and N oc t^^,., to the much steeper, GW-driven N oc tvLv - below 
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Figure 4. Left Panel: The figure shows the time Tqw that a SMBHB, whose 
coalescence is driven entirely by gravitational radiation, spends at a radius where 
the orbital time is within torb- Both Tqw a^nd to^b are in redshifted units (as 
measured on Earth). Each of the shaded regions corresponds to a different total 
BH mass ([1 + z]Mtot = [10^'^' ''■*'^'l°]M0; top to bottom), and shows three 
different mass ratios (q = [1, O.f , 0.01]; bottom to top, or solid, short-dashed, and 
long-dashed, respectively). The quantity on the y axis can be interpreted as the 
duty— cycle for each source during which it exhibits periodic variability on the 
time-scale tvar ~ toib- The thick (red) portion of the lines and the vertically 
shaded (red) areas indicate the orbital separations, for each mass, where binary 
evolution may no longer be dominated by GWs. The slope of the arrows indicate 
the fiatter evolution of the binary when it decays due to viscous torques: the 
flatter arrow corresponds to the slope T oc t^^^^^ for near-equal mass binaries, 

and the slightly steeper arrow to T oc t^j^^^ for ij <SC 1 (see Appendix in ref. [10] 
for details). For reference, the x axis labels on the top show the orbital radius 
corresponding to each orbital time; for M = IO^Mq (or M7 = 1) and z = 0, the 
plotted range covers 1.6 X lO-"^ - 7.5 X 10-=* pc, or 17-77,000 Schwarzschild radii. 



Right Panel: This panel shows the sky coverage of a future survey, required to 
find at least one periodic source, with a period of at most 20 weeks, as a function 
of the i-band magnitude limit on the variable component. This will ensure that a 
multi-year survey, could find many more sources with longer periods, e.g. periods 
up to 1 year, demonstrating the periodic nature of the variations, and allowing 
to deduce the abundance of periodic sources as a function of period. The three 
solid curves show results assuming that the variable component is a fraction 0.3, 
0.03, and 0.003 of the Eddington luminosity of the binary. Each dashed curve 
shows the corresponding mass of the SMBHB of which there would be a single 
example of a fvar = 20-week periodic variable in the survey (~ 20 sources with 
the same luminosity are then predicted with a tvar = 60 week period in the same 
survey volume, assuming pure GW— driven evolution, or fewer such longer-period 
sources, if viscosity speeds up the evolution significantly). 

a characteristic ivar ^ tens of weeks (corresponding to ^ 10^ Schwarzschild radii for 
equal-mass ~ IO^IVIq binaries), could confirm (1) that the orbital decay for sources 
below a characteristic ivar is indeed driven by GWs and (ii) also that circumbinary 
gas is present at small orbital radii and is being perturbed by the BHs - serving as a 
proof of concept for eventually finding actual LISA counterparts. 
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Finally, we mention two other possibilities to identify coalescing SMBHBs before 
LISA'S launch. First, [24] and [25] followed the response of the gas disk around a 
recoiling SMBHB, similar to our calculation described in § 3, but on much longer 
timescales (~ 10^ years). Ref. [24] found that the shocked gas, thrown out of the 
disk plane by an oblique kick, may produce bright flares in X-ray lines (assuming it 
remains optically thin). Ref. [25], on the other hand, assumed that the disk gas is 
optically thick, and that the shocks are promptly dissipated in the disk plane, and 
derived a characteristic infrared light curve. The above effects arise from the change 
in the gravitational potential, following the recoil, i.e. from the "shaking" of the 
disk. A different possibility is that once the BHs have merged, the torques from 
the rotating quadrupolar potential cease to act on the gas outside the inner cavity, 
allowing the gas to accrete onto the merged binary on the viscous time-scale [23] . An 
X-ray "afterglow" may then occur after a delay of ~7(1 -|- z){M /W^Mq)'^-^'^ years. 
This could be detectable in LISA follow-up observations, or perhaps without a LISA 
trigger for very massive BHs at high z, where the time-scale is long. 

Acknowledgments 

ZH thanks the organizers of the conference. The work described here was supported by 
NASA (grant NNX08A1135G), and by the Polanyi Program of the Hungarian National 
Office for Research and Technology (NKTH). 

References 

[I] P. J. Armitage, & P. Natarajan, Astrophys. J. 567, L9 (2002). 
[2] P. Artymowicz, & S. H. Lubow, Astrophys. J. 421, 651 (1994). 

[3] P. Artymowicz, & S. H. Lubow, Astrophys. J. Lett. 467, L77 (1996). 

[4] J. E. Barnes, & L. Hernquist, Ann. Rev. Ast. & Astrophys. 30, 705 (1992). 

[5] M. Campanelli, et al., Phys. Rev. Lett. 98, 231102 (2007). 

[6] L. Corrales, A. MacFadyen, & Z. Haiman, Astrophys. J., to be submitted (2008). 

[7] J. Cuadra, P. J. Armitage, R. D. Alexander, & M. C. Begelman, Mon. Not. Royal Ast. Soc., 

submitted; arXiv:0809.0311 
[8] C. Deflfayct, & K. Mcnou, Astrophys. J. 668, L143 (2007). 

[9] A. Escala, R. B. Larson, P. S. Coppi, & D. Mardoncs, Astrophys. J. 607, 765 (2004). 
[10] Z. Haiman, B. Kocsis, & K. Mcnou, Astrophys. J., submitted; arXiv:0807.4697 

[II] D. E. Holz, & S. A. Huglics, Astrophys. J. 629, 15 (2005). 

[12] P. B. Ivanov, J. C. B. Papaloizou, & A. G. Polnarcv, Mon. Not. Royal Ast. Soc. 307, 79 (1999). 

[13] B. Kocsis, Z. Frci, Z. Haiman, & K. Mcnou, Astrophys. .]. 637, 27 (2006). 

[14] B. Kocsis, Z. Haiman, K. Mcnou, & Z. Frci, Phys. Rev. D 76, 022003 (2007). 

[15] B. Kocsis, Z. Haiman, & K. Mcnou, Astrophys. J. 684, 870 (2008). 

[16] B. Kocsis, & A. Locb, Phys. Rev. Lett. 101, 041101 (2008). 

[17] R. N. Lang, & S. A. Hughes, Phys. Rev. D 74, 122001 (2006). 

[18] R. N. Lang, & S. A. Hughes, Astrophys. J. 677, 1184 (2008). 

[19] Z. Lippai, Z. Frci, & Z. Haiman, Astrophys. J. Lett. 676, L5 (2008). 

[20] A. MacFadyen, & M. Milosavljevic, Astrophys. J. 672, 83 (2008). 

[21] P. Madau, "Cosmic Star Formation History and the Brightness of the Night Sky", in Prom 
Extrasolar Planets to Cosmology: The VLT Opening Symposium, edited by J. Bergeron and 
A. Renzini (Springer- Verlag, 2000), p. 52.; arXiv:astro-ph/9907268 

[22] F. Melia, & H. Falcke, Ann. Rev. Ast. & Astrophys. 39, 309 (2001). 

[23] M. Milosavljevic, & E. S. Phinncy, Astrophys. J. 622, L93 (2005). 

[24] G. A. Shields, & E. W. Bonning Astrophys. J. 682, 758 (2008). 

[25] J. D. Schnittman, & J. H. Krolik, Astrophys. J. 684, 835 (2008). 

[26] B. F. Schutz, Nature 323, 310 (1986). 

[27] J. H. Taylor, & J. M. Weisberg, Astrophys. J. 253, 908 (1992). 
[28] A. Vecchio, Phys. Rev. D 70, 042001 (2004). 



